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ABSTRACT 


Localized velocity, temperature, and species concentration measurements in rocket flow fields are needed to evaluate 
predictive computational fluid dynamics (CFD) codes and identify causes of poor rocket performance. Velocity, 
temperature, and total number density information have been successfully extracted from spectrally resolved Rayleigh 
scattering in the plume of small hydrogen/oxygen rockets. Light from a narrow band laser is scattered from the moving 
molecules with a Doppler shifted frequency. Two components of the velocity can be extracted by observing the scattered 
light from two directions. Thermal broadening of the scattered light provides a measure of the temperature, while the 
integrated scattering intensity is proportional to the number density. Spontaneous Raman scattering has been used to 
measure temperature and species concentration in similar plumes. Light from a dye laser is scattered by molecules in the 
rocket plume. Raman spectra scattered from major species are resolved by observing the inelastically scattered light with a 
linear array mounted to a spectrometer. Temperature and oxygen concentrations have been extracted by fitting a model 
function to the measured Raman spectrum. Results of measurements on small rockets mounted inside a high altitude 
chamber using both diagnostic techniques are reported. 

2 INTRODUCTION 


Improving the performance of on-board rockets yields rewards in terms of extended mission life, or enhanced mission 
goals through a reduced propellant load. 1 Computational fluid dynamics (CFD) codes are routinely used both to evaluate 
rockets which are under development and to establish the effect of design changes on existing rockets. In some cases an 
acceptable prediction of rocket performance can be made with the current level of predictive technology. For high and 
medium thrust rockets, predictions of the global parameters such as thrust and specific impulse, show good agreement with 
measurements. Using these codes for predicting the performance of certain low thrust rockets, however, leads to inaccurate 
predictions, the specific causes of which are not identified. 2,3 Improved predictions for these rockets can only be made based 
on a better understanding of the physics of small rockets, a refined analysis of the individual rocket, and an increase in 
numerical capabilities. Such developments are aided by the availability of localized fluid dynamic and thermodynamic data 
upon which these improved codes can be verified. 

The rockets tested in this investigation were gaseous hydrogen/oxygen rockets mounted inside a high altitude simulation 
chamber. The absence of ambient light inside the high altitude chamber combined with the low luminosity of the water 
vapor rich plume provided ideal conditions for spectroscopic probing by lasers. Laser Doppler anemometry requires the use 
of particles. The large gradients inherent in small rockets as well as the difficulties in introducing particles precluded the 
application of this technique to small rockets. This prompted the development of a technique based on Rayleigh scattering. 
This technique yields velocity and in some cases the temperature and number density of the molecules in a high velocity 
plume. 4 For Rayleigh scattering a narrow band laser beam (pulsed or cw) was directed through the plume. The laser light 
scattered from the moving molecules was Doppler shifted and collected at two locations. The magnitude of the observed 
Doppler shift is dependent on the relative angles between the incident and scattered light and the gas velocity. Since the 
relative angle between the incident and collection optics was known, the mean velocity of the observed gas volume was 
determined from the magnitude of the Doppler shifts measured at each location. 



In addition to a mean gas velocity, each molecule possesses its own thermal motion which depends on the gas 
temperature and the molecular mass. When the gas is in thermal equilibrium, the molecular velocity obeys a Maxwellian 
distribution, which results in a Gaussian profile for the scattered laser light. 5 Therefore, in the absence of other significant 
broadening parameters, the degree of broadening for a single species can directly be related to the temperature 
Accordingly, by resolving the line shape of the scattered light, and relating this to the reference laser line shape, it was 
possible to directly estimate the molecular velocity and translational temperature. Furthermore, since the scattering intensity 
for a single species is linearly proportional to the molecular number density, total number density could be estimated for the 
mixture by comparing the total scattered intensity and the intensity scattered by a reference gas. 

Spontaneous Raman scattering was implemented for measuring temperature and species number densities. Light from a 
pulsed dye laser was focused in the rocket plume. Both Rayleigh and Raman scattering processes occur inside this focal 
volume. Raman scattering occurs when light photons exchange energy with rotational and/or vibrational internal molecular 
energy modes. Incident photons that lose energy to a molecular energy system are shifted to a lower spectral frequency. This 
spectral shift depends on the scattering species since these energy losses are discrete and equal to a change in molecular 
energy quantum. Thus each observed shifted spectral line indicates the presence of a specific species, characterized by the 
magnitude of the shift from the incident laser line. For a set of molecules in thermal equilibrium, the population over all 
possible energy levels can be described by a Boltzman distribution. Since quantum jumps are equally likely from all initial 
energy levels, the distribution of the Raman scattered photons over the different possible shifted lines for each species 
reflects this Boltzman distribution. Light scattered from the probe volume was collected and spectrally separated 
Measurement of the intensity distribution of this light over the possible spectral range for the measured species yielded the 
Boltzman distribution from which the temperature was determined. The overall intensity of this distribution depends 
linearly on the incident laser intensity and on the number density of the scattering species. Because the laser intensity was 
known, the number density was directly determined from the scattered intensity. 

Both techniques are impeded by the presence of ambient light. Raman scattering is a naturally weak process which 
requires minimum ambient and stray light to achieve acceptable signal to noise ratio (SNR). Rayleigh scattering is three 
orders of magnitude stronger and is therefore less sensitive to experimental conditions. However, for the measurements 
reported in this paper, the combined effects of low scattering number densities and the need for highly resolved spectra 
resulted in low signal levels. Judicious placement of pinholes, baffles and beam stops was therefore required to achieve an 
acceptable SNR. The limited number of windows of the high altitude simulation chamber were blocked to remove all 
ambient light, and extensive use was made of optical fibers for accessing the chamber. 

Velocity and temperature measurements based on spectrally resolved Rayleigh scattering were successfully made on a 
lab model gaseous hydrogen/oxygen rocket engine with a nozzle area ratio of 30:1. The nozzle was modified for the Raman 
measurements since the number densities in the exit plane of the full scale rocket were prohibitively low. A reduced area 
ratio nozzle was installed which made number density and temperature measurements feasible in the exit plane. The results 
obtained with both techniques were compared to predictions by existing CFD codes. Results are reported in this paper. 

3, ROCKET TEST FACILITY 


The rockets tested were designed for space operation. To simulate a space environment, tests were conducted in a high 
altitude simulation chamber.** This chamber is 0.91 m in diameter and 1.81 m in length. The chamber is evacuated by 
means of a supersonic diffuser which is driven by a two stage venturi ejector stack. The 0.102 m diameter diffuser fully 
captured the rocket plume The diffuser performance was optimized by locating the diffuser entrance from 0 to 0.10 m 
behind the rocket exit plane When optical access was required, a trade-off between the diffuser distance and optical access 
was necessary. Under optimal conditions, chamber pressure was maintained at 1.3 kPa. simulating conditions at a height of 
35 km. During tests when the plume was not captured sufficiently, the pressure gradually rose. Tests were usually aborted 
after the pressure exceeded 6 9 kPa. A detailed description of the high altitude simulation facility and the optical access is 
given in Ref. 6. 

The rockets, shown schematically in Fig. 1 were in the 1 10 N thrust class. 2 They were designed to operate with gaseous 
hydrogen and oxygen with an oxygcn-to-fucl ratio of 3 .0-5.0. The propellant was supplied by high pressure trailers. A 
combination of upstream regulators and calibrated sonic venturis controlled the propellant flow rates to the thruster which 
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were continuously monitored. The 0.025 m diameter, 0.102 m long combustion chamber converged into a nozzle with a 
throat diameter of 0.0128 m and exit area ratio of 33:1. Gaseous hydrogen was used to regeneratively cool the combustion 
chamber and nozzle. The hydrogen entered wall passages through a channel in the nozzle exit plane and propagated 
upstream through the chamber walls toward the combustion chamber. The hydrogen flow was then split by means of 
calibrated flow passages. Part of the hydrogen was directed towards the core flow of the injector and part was diverted to a 
sleeve which distributed the fuel in a thin film over the combustion chamber wall. The amount of fuel diverted for fuel film 
cooling was regulated between 50% and 75%. The thruster injector was a stacked platelet injector. Oxygen was injected 
upstream of a spark plug and flowed past the spark plug tip towards the combustion chamber At the spark plug tip, 
hydrogen was injected through radial passages in the stacked platelets into the oxygen flow and was entrained in the 
recirculation zone created behind the blunt spark plug tip. After initially igniting the oxygen rich mixture, the spark plug 
functioned as a flame holder, creating turbulent mixing and maintaining combustion. The oxygen rich product mixture 
entered the combustion chamber and gradually consumed the hydrogen fuel film at the wall. 


Several conditions were favorable for applying Rayleigh scattering to the plume. The choice of propellants led to a 
plume with extremely low luminosity in the wavelength range of interest for Rayleigh scattering, providing an acceptable 
signal-to-noise ratio. Furthermore, the low' pressure in the exit plane guaranteed that the mean-free-path between molecules 
was much larger than the scattering interaction wavelength. Accordingly it could be assumed that the scattering occurred in 
the "collisionless" regime, where uncorrelated scattering of molecules could be assumed. This greatly simplified the 
scattering analysis. 7 

The Rayleigh scattering analysis assumed water to be the single scattering species. Because the overall mixture ratio was 
fuel rich it is reasonable to assume that the plume was fuel rich. Water molecules have a much larger Rayleigh scattering 
cross section than hydrogen and the much lower molecular mass of hydrogen causes its scattering contribution to be spread 
over a much wider spectral range. Therefore water vapor was assumed to be the dominant scattering species. Due to the fuel 
film cooling however, there was a considerable variation in local mixture ratio, gradually varying from an oxygen rich core 
to a fuel rich perimeter region. As previously stated, water vapor was the dominant scatterer in those regions that were fuel 
rich. Oxygen contributed significantly to the total scattering for measurements in the core. This made it impossible to 
estimate the temperature and total number density. 


The extremely low luminosity of the hydrogen-oxygen combustion products in the plume was also favorable for Raman 
scattering. As with the Rayleigh scattering, the low density made a simple, collision free analysis possible. However, the 
same low density made it difficult to detect Raman scattering in the exit plane during the 30 second test times. Therefore, 
for Raman scattering experiments, the nozzle length was reduced to an area ratio of 1.86:1. A simple one-dimensional 
nozzle analysis showed the pressure in this exit plane to be about 0.1 MPa., theoretically providing a sufficient number 
density for Raman scattering diagnostics. The redesigned test rocket wall was water cooled. The injected hydrogen which 
provided film cooling was lower in temperature than the hydrogen injected in the regeneratively cooled rocket. As a result, 
the fuel film injection velocity was different than during the Rayleigh experiments. 
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Fig. 1. Regeneratively cooled hydrogen-oxygen thruster flow field with fuel film cooling. 
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4. SPECTRALLY RESOLVED RAYLEIGH SCATTERING 


Rayleigh scattering occurs when an incident electromagnetic wave and a molecule interact without exchanging energy 
Three conditions have to be satisfied for a relatively simple treatment of both theory and application. The wavelength of the 
incident wave should be much larger than the size of the scattering molecule. If this is satisfied, a simple oscillating dipole 
is induced inside the molecule, which creates its own electromagnetic field and subsequent radiation (scattering). The dipole 
oscillates with the frequency of the incident wave, and consequently the scattered wave is of the same frequency For a 
simple "collisionless" treatment of Rayleigh scattering theory, the scattering number densities have to be low enough to 
ensure that the collisional mean free path of the molecules is much larger then the interaction wavelength If that is the 
case, the incident wave will interact with the molecules in an uncorrelated manner 7 and the resulting Ray leigh scattering is 
characterized by a simple Gaussian shape. Finally, the scattering volume should consist of a single component gas or the 
mixture composition should be known. When a wave is incident on a single component gas volume, the scattered intensity 
depends only on the molecular number density, and can be used for measuring number densities. Since the strength of the 
dipole that is created varies with type of molecule, the scattered intensity varies with type and number of scattering 
molecules. In addition, because the bandwidth of the Gaussian shape depends on the molecular weight, the contribution or 
each species to the final profile is a Gaussian of widely varying shape. It is therefore generally difficult to measure 
molecular number densities in a gas mixture if the gas composition is not known a priori. 

The directionality of Ray leigh scattering depends on the polarization of the incident wave. For a linearly polarized 
incident wave, the dipole oscillates one-dimensionally and the radiation shows a toroidal distribution, with the strongest 
scattering perpendicular to the oscillating dipole. It is therefore important to detect the Rayleigh scattered radiation 
perpendicular to the electric field vector. 


4, 1 Theory 

The random velocities of the individual molecules of a thermalized gas volume at rest follow a Maxwellian velocity 
distribution. Because of this, the Rayleigh spectrum that results when an incident plane wave of single frequency is 
scattered by this gas volume exhibits a Gaussian shape centered around the incident wavelength. This shape possesses a 
bandwidth (FWHM) which is directly proportional to the square root of the gas temperature. For a gas volume with mean 
velocity not equal to zero, the peak position of the observed spectrum is shifted from the incident wavelength. The shift 
depends on the gas velocity and the relative orientation between the incident and observed light, as illustrated in Fig. 2a for 
a backscattcr and Fig. 2b for a forward scatter geometry-. A detailed derivation shows that the power scattered by a single 

component gas in the direction of observation given by "scattered" vector k s into solid angle dO in frequency interval df 


P s (f)dfdQ = IonV^ 


dcr 

dD 


sin 2 ^ S(f) dfdQ 


0 ) 


where I is the incident irradiancc in W/m 2 . n the molecular number density in nr 3 , and V sc is the scattering volume The 


differential Rayleigh scattering cross section 


dcr 

dQ 


is given in m 2 /sr, and X ' s th e angle between the incident electric field 


vector k„ and k s The frequency distribution function S(f) is the normalized spectrum of Rayleigh scattered light. This 
normalized spectrum can be obtained by integrating the Maxwellian velocity distribution over the velocity space. This 
gives: 4 


S(f) df = 


2 v/*T ^-(2,rf - K ii) 2 /a 2 K 2 
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where a = ( 2K'T/m) 12 is the "most probable speed" of the molecules in m/sec, and K is the scattering vector 

magnitude, where the scattering wave vector can be obtained from the relation K “ k s - k 0 . The thermodynamic 
constant K is Boltzmann’s constant (1.38 x 10‘ 23 J/K), T is the gas temperature in Kelvin, and m is the molecular mass in 
kg. The vector D indicates the mean velocity whose magnitude is given in m/sec. Eq. 2 represents a Gaussian vclocitv 
distribution with a bandwidth (FWHM) of 0.265 a K Hz. 



Fig. 2. Vector diagram of Rayleigh scattering; (a) backscatter and (b) forcscatter. 

For a single component gas with known molecular mass, the temperature can be extracted from the bandwidth of the 
observed scattered light. The mean velocity of the molecules can be found from the location of the peak of the Gaussian 
profile with respect to the incident wavelength as is shown in Fig. 3. The number density can be estimated from the total 
integrated Rayleigh scattered light. To obtain a density measurement, a reference measurement is necessary, where the total 
integrated Rayleigh scattering is measured with identical optical geometry on a gas with known number density and know n 
differential scattering cross section. The number density measured is linearly proportional to the ratio of the total integrated 
Rayleigh scattering powers for both cases. The absolute number density can be found from the reference number density and 
ratio of differential scattering cross sections. 

For a multi-component system with unknown composition, the total scattered Rayleigh profile can be assumed to be 
made up of the sum of the contributions from each individual species. If all species possess the same mean velocity', each 
contribution consist of a Gaussian lineshape located at the same spectral location with different spectral widths, dependent 
on the scaltcrers' molecular mass. The determination of temperature and number density from such a profile requires 
additional information However, in many cases, a single species dominates the scattering process and a simple single 
component curve fit can yield the desired data. 

4 2 Application. 

Plumes of hydrogen -oxygen rockets consist predominantly of water vapor. Both pulsed and continuous wave (cw) 
strategies were employed to generate Rayleigh signals. The cw p technique utilized the 514.5 nm line of a cw Argon-Ion 
laser. The pulsed strategy utilized the second harmonic output from an injection seeded Nd:Yag laser (X = 532 nm, 150 
MHz bandwidth). Both approaches satisfied the condition that the wavelength of the incident wave was much larger than 
the 1 nm effective diameter of the watermolecule. The light was vertically polarized and the detection direction, given in 

Eq. I as the vector k s , was in a horizontal plane. This implies that the angle % was 90°, which maximized the scattered 
power in the detection direction The optical geometry for the cw strategy is shown in Fig. 4. Various steering optics 
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directed the laser beam through the optical port at 60°, parallel to the optical table mounted underneath the thruster, and 
between the diffuser and thruster. A 1000 mm focal length lens was used to focus the beam to a diameter of 180 pm in the 
near-field region of the thruster. Since stray laser light could impede the measurement by entering the detection system, 
extensive use was made of baffles, masks, high quality' optical windows and beam stops. For the pulsed strategy, the laser 



Fig. 3. Typical spectral scan of Rayleigh scattering over two times the free spectral range. 


Two fiber optic collection probes were positioned at angles 8 S = 31° and 149° (see Fig. 2) for fore- and backscatter, 
respectively, to determine two components of the velocity. Each probe consisted of two lenses, one lens to collect scattered 
light and a second one to refocus it into a fiber, and a 10 m long fused silica fiber with a diameter of 1000 pm which gave 
an effective collection f# of 4.55. This geometry resulted in an adequate collection solid angle without introducing aperture 
broadening effects caused by the variation of collection angle over the collection optics. The fiber optic guided the scattered 
light through the altitude chamber bulkhead to an optical table, where the fiber output was collimated with an achromatic 
doublet and passed through a single pass scanning Fabry-Perot interferometer. The free spectral range of this interferometer 
was 21.4 GHz, sufficient to resolve the bandwidth (FWHM) of a nominal Doppler broadened Rayleigh spectrum of water 
(T«1000 K. X=514.5 nm, and 0 S =31°), which was approximately 6 GHz. The filtered light leaving the Fabry-Perot was 
focused through a 1000 pm aperture; collimated; passed through a multilayer dielectric interference filter to eliminate 
luminous background contributions; and detected with a 20% quantum efficiency, Gallium-Arsenide photomultiplier tube 
(PMT). A single threshold photon counter was used to measure the resulting Rayleigh scattering signal. Piezoelectric spacer 
elements were used to vary' the Fabry-Perot mirror spacing over two adjacent orders of interference, thereby tuning the 
delected PMT signal through the Rayleigh spectrum. A typical spectral scan is displayed in Fig. 3. 

Temperature and velocity were extracted from the measured Rayleigh spectrum by curve fitting a model function based 
on Eq. 1 to this spectrum. It was necessary to expand the model function to include the instrument response function. This 
instrument function was experimentally determined for both the unshifted, unbroadened argon-ion laser line, and the 
Nd:Yag second harmonic output. A small portion of laser light, reflected from the altitude chamber window was collected, 
attenuated with neutral density filters and coupled into a reference fiber. This fiber guided the light into the test section 
where it was mixed with the Rayleigh scattered light and collected with the collection probes. This reference signal could be 
turned on and off with the placement of a shutter on the collection side of the reference fiber Excellent agreement was 
obtained between the experimentally measured and calculated reference spectra. For number density measurements, 
calibration spectra were obtained by traversing the measurement probe volume across the rocket exit plane while the facility 
was purged with a nitrogen gas of known temperature, density, and scattering cross section. High temperature spectra in the 
rocket plume were obtained at the same radial locations by traversing the probe volume along the path of the laser beam 
during a rocket fire. 
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Use of finite diameter optical fibers reduced the spectral resolution of the interferometer due to the increase in solid 
angle from the light collimated through the instrument. 7 This limited the sensitivity of the temperature measurement. Eq 1 
shows that the Doppler width was proportional to the magnitude of the scattering vector This vector was much larger in the 
backscattcr geometry, which caused the Doppler broadening to be much greater in the backscattcr direction, surpassing the 
sensitivity of the instrument even at low temperatures (a theoretical limit of 245 K). As a result, temperature measurements 
could only be made in the backscatter geometry. 



Fig. 4. Diagram of Rayleigh scattering optical arrangement with cw laser source. 

An error analysis based on bias errors and the maximum likelihood parameter estimation procedure used in the curve 
fitting process, showed that the typical uncertainty in the mean velocity estimate was 200-300 m/s, independent of the 
magnitude of the mean velocity. Temperature determination was slightly less accurate, primarily due to the sensitivity of the 
estimation to noisy spectral profiles. Typical temperature uncertainty was 100-150 K. Number density estimates for a single 
component gas typically fall within 2-5% of the mean value. 

5. VIBRATIONAL-ROTATIONAL RAMAN SCATTERING 

Raman scattering is a process that occurs simultaneously with Rayleigh scattering when a plane electromagnetic wave is 
incident on a molecule with vibrational and rotational degrees of freedom. Because the internal molecular energy field 
changes during molecular vibration, the dipole induced by an incident wave changes accordingly. If the frequency of the 
incident wave is much higher than the molecular vibrational frequency, this change only becomes perceptible across many 
oscillations of the incident electric field. As a consequence, Rayleigh scattering which results from a constant oscillating 
dipole will be the dominant feature of the scattering spectrum. But a modulation of the dipole occurs with a period equal to 
the vibrational period of the molecule which results in the appearance of a much weaker spectral line shifted away from the 
Rayleigh frequency by a magnitude equal to the vibrational frequency. A similar spectral line appears due to the fluctuation 
in induced dipole during a molecular rotation. In this case, the period of change is half the period of one rotation since the 
molecular susceptibility towards dipole generation is equal for a molecule rotated 180°. The lines are therefore shifted away 
from the incident frequency by twice the rotational frequency. 

A quantum mechanical description of this phenomenon is based on the exchange of energy between an incident photon 
and rotational-vibrational molecular energy modes. Since a molecule can only exist in discrete vibrational and rotational 
energy states, it absorbs enough of the energy of an incident photon to cause a transition to a upper energy state. The 
resulting ’'scattered” photon is less energetic and appears in the spectrum shifted to a lower frequency. Or in the case that 
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the molecule is already in an upper state, it might release enough energy to the photon to return to a lower energy state. It 
then scatters a more energetic photon which appears in the spectrum shifted to a higher frequency. Wave function 
calculations 8 show that vibrational transitions are only allowed between adjacent levels, which means a change of one 
vibrational quantum number A vibrational transition to a upper energy level, given by a quantum change of Av = +1, is 
referred to as Stokes scattering. A vibrational transition to a lower level, indicated by quantum change Av = -1, is referred 
to as anti-Stokes scattering. Similar calculations show that rotational transitions are only allowed between every other level, 
thus corresponding to a rotational quantum number change of two. Scattering caused by rotational transitions given by 
quantum jump AJ = +2 are referred to as S-branch scattering. A change of AJ = -2 yields the O-branch and the pure 
vibrational transition AJ = 0 gives the Q-branch. 

The energy of the scattered photon, and accordingly the spectral location at which it appears, depends on the amount of 
energy that has been exchanged with the molecule. This is equal to the difference between the two energy levels. Since these 
differences are unique for each molecule, the spectral line location can be used to identify the specific species that cause the 
scattering. 

5. 1 Theory 

The distribution of the molecules over all possible energy levels of a gas volume in thermal equilibrium is described by 
the Boltzmann distribution. Since each of these molecules contributes to the scattered light, many lines appear in the 
spectrum that are shifted from the incident light. Vibrational and rotational motions of the molecules interact, with the 
magnitude of interaction dependent on the specific levels. Thus the amount of energy that is needed for a pure vibrational 
transition (Q-branch) of a single molecule, depends on the initial rotational level of that molecule. As a consequence, each 
rotational level in a pure vibrational transition appears at a different wavelength in the spectrum. The strength of each of 
these lines is directly proportional to the number of molecules in that specific rotational energy mode. The resulting 
distribution reflects the initial rotational distribution, which is an accurate indication of the rotational temperature. Because 
the differences in line locations are small, high resolution is needed to resolve the distribution. As the result of poor 
resolution and line broadening, they usually show up in the spectrum as one single line. Spectral lines caused by 
simultaneous rotational-vibrational transitions (O- and S-branches) lie further apart. However, these transitions have a 
lower transition probability, given by the Placzek-Teller coefficients. Accordingly, they are easier to resolve but more 
difficult to detect. 


At higher temperatures, more vibrational levels become populated. As a result, upper vibrational transitions (v=l— >2, 
v=2-»3) cause new series of lines to appear, each representing the full rotational distribution of the particular vibrational 
state. The relative strength of each v ibrational series reflects the molecular number density in that initial vibrational energy 
mode. Comparison of the different vibrational series line intensities gives the vibrational temperature. Comparing the total 
integrated scattering intensity over all transitions to that of a reference gas at known number density and temperature yields 
the species number density. 

For a single component gas in thermal equilibrium, the total power scattered in the direction of observation in solid 
angle dQ as the result of one allow able transition can be given as: 9 


dn 


do = I 0 n V_ 

J0. v\j'->v # .r 

where I 0 is Ihc incident irradiancc in W/m 2 . fl the molecular number density in nr 3 , and is the scattering volume in 


dcr 
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m 3 . The differential Raman scattering cross section 
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is given in 


m 2 /sr. 
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The shape function S v J , reflects 


the Boltzmann distribution of the molecular population in its initial undisturbed state. This function only depends on the 
temperature, and is given by: 
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where g n is the nuclear spin factor. This factor is a weighting factor which reflects the naturally occurring difference in 
spin populations between the even and odd rotational distributions (ortho- and para- states) of diatomic homonuclear 
molecules and is 1 for heteronuclear molecules. For example, for hydrogen, the even:odd ratio is 1 :3, for nitrogen, this ratio 

is 2:1, and for oxygen, the even states are missing and its nuclear spin factor g n = 0 for the even lines, which means that 
transitions starting from even rotational energy levels do not exist. The partition function Z is the rotational -vibrational 

probability distribution of the molecules over all possible energy states and the energy E v , y is the rotational-vibrational 

energy of the initial molecular state. A single prime, \ refers to the initial molecular state and the double prime, ", to the 
final state. 


The differential Raman scattering cross section depends on the type of scattering molecule, the specific rotational- 
vibrational transition, and the fourth power of the frequency of the scattered wave. It also depends on the angle 0 between 
the incident and the scattered radiation, on the polarization of the incident and observed radiation, and on the angle 
between the incident electric field vector and the direction of observation. For any specific experiment, the last three species 
independent properties remain constant. The dependence of the scattering cross section for a particular species is on the 
wavelength of scattering and the specific transition Measured values of the scattering cross sections for different species are 
usually given for a specific optical configuration and incident wavelength. They are tabulated as the ratio between this cross 
section and the scattering cross section of the nitrogen Q-branch, 10 for which the absolute values are well documented 
Corrections for different experimental optical characteristics and wavelength can be implemented accordingly. 11 12 

It is often not necessary to know the absolute scattering cross section as long as the different dependencies arc accounted 
for. The type of species can be determined from the spectral location where scattered light is observed, and the temperature 
can be determined from the intensity distribution of the scattered spectrum. The underlying profile of this distribution 
reflects the Boltzmann population and as such is a direct thermometer. To extract an accurate value for the temperature, the 
shape of the basic distribution has to be corrected for its dependence on the rotational-vibrational quantum numbers. The 
change in scattering cross section as a function of transition type and scattering wavelength has to be accounted for The 
functional dependence is: 


00 P J.AJ ( V ' + l ) K - ^Ks) 4 ( 5 ) 

o.v\y-*yf,r 

where Pj AJ are the rotational quantum number dependent Placzek-Tellcr coefficients. 11 A different J-dcpcndcnt coefficient 
exists for each of the rotational branches (0-,Q- and S-branch). The term (v f + I) indicates the vibrational dependence of 
the line strength of the Stokes transition. This term becomes (v*) for the anti-Stokes transition, indicating that from the 
vibrational ground state (v ,= = 0) such a transition does not exist. The term v 0 is the frequency of the incident wave and 

v RS is the Raman shift. This shift is the energy difference between the initial and final energy state of the transition divided 
by Planck’s constant, and is species dependent. 

After the temperature has been extracted, the species number density can be determined. The ratio of the total scattered 
power versus the total scattered power of the nitrogen Q-branch obtained under known conditions gives the number density 
To obtain high accuracy, the number density and temperature of a reference sample of the measured species must be 
carefully controlled. This sample can then be used as a secondary calibration standard 


dcr 

dQ 
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5 2 Application 


A full description of the Raman diagnostics facility is given in Ref. 14. The layout of the facility is displayed in Fig. 5 
A flashlamp pumped dye laser, lasing at a nominal wavelength of 595 nm, was directed across the plume to generate 
spontaneous Raman scattering. The pulse length was approximately 2 psec and the nominal pulse energy about 1 
Joule/pulse. The bandwidth (FWHM) of the incident radiation was 0.03 nm (or 0.9 cm 1 ). A pulsed input was chosen to 
minimize the background luminosity by gating the detection equipment. Use of a shorter pulse length and/or shorter 
wavelength incident radiation created concern about laser induced gas breakdown.* 5 For wavelengths in the ultraviolet 
region, fluorescence signals exceeded the Raman scattering signals. Furthermore, at longer wavelengths, individual Raman 
transition lines arc further apart, allowing a better resolution of the band profile. An additional advantage of a wavelength 
in the visible is that conventional optical components could be used for transmission and manipulation. However, because 
the scattered intensity scales with the fourth power of the scattering wavelength the signal levels are smaller at longer 
wavelength ranges. The decay in laser pulse energy' as a function of time was continuously monitored by splitting 1% of the 
pulse energy off and directing it to a pulse energy meter 

The pulsed beam was coupled into an 800 pm optical fiber of 20 m length which guided it through the high altitude 
facility bulkhead. The pulse emanating from the fiber was collimated and refocused in the exit plane of the rocket nozzle 
The optical axis lay in a plane perpendicular to the thruster axis. Receiving optics with an f# of 2.4 were mounted in the 
same plane perpendicular to the incident optics. These optics consisted of a collection lens and a refocusing lens w hich 
coupled the detected scattering into a 20 m length, 800 pm diameter optical fiber. This fiber guided the scattered light 
through the altitude chamber bulkhead into a diagnostics room. The scattered light emanating from this fiber was 
collimated and filtered with a Raman notch filter with a better than 10 6 signal attenuation at the laser line. The light was 
focused into a 0.5 m spectrometer with 300 or 1200 groove grating, dependent on the application. For simultaneous 
multispecies detection, the 300 groove grating gave a 6.4 nm/mm spectral resolution. This allowed simultaneous detection 
of oxygen, hydrogen, nitrogen and water with a single 25 mm array detector. For temperature measurements where line 



Two detection schemes were used. A 700 pixel linear array with 25pm pixel size was gated such that the gate was 
coincident w ith the laser pulse. With the use of background subtraction, this pulsed system could extract very weak signals 
from the noise. The measured temperature could be extracted from the measured line shapes using a model function based 
on Raman theory and the instrument function of the detection equipment. Calibration factors were obtained from a 
reference gas at known temperature and number density. A maximum likelihood parameter estimation procedure completed 
the data reduction The instalment function was measured using a helium neon laser line and a mercury lamp. Oxygen and 
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nitrogen spectra obtained in atmospheric air were used as calibration spectra. Typical nitrogen spectra, both measured and 
calculated, for two different temperatures are shown in Fig. 6. The spectrometer entrance slit width for these measurements 
was 200 pm. 

At very low number densities and limited test time, the linear array did not provide sufficient sensitivity. Therefore a 
GaAs photomultiplier tube with a higher quantum efficiency (Q.E.=12% @ 700 nm) was used. A boxcar averager and 
integrator was used to gate and integrate the PM tube signal with optimum noise rejection as is shown in Fig. 5. With an 
exit slit of 200 pm a nominal spectral resolution of 0.33 run was obtained. By scanning the spectral range of interest and 
integrating at each spectral location over 300 laser pulses, a spectral profile could be obtained. For very weak signals it was 
necessary to open the spectrometer entrance slit further to allow more light. The 800 pm collection fiber had a Numerical 
Aperture (N.A.) of 0.22. Collimating and refocusing optics projected a 800 pm image of the collection fiber exit on the 
spectrometer entrance slit, slightly over filling the (110x110 mm) grating. For maximum light transmission, the entrance 
and exit slit were set to 800 pm. The resolution of this method however was significantly worse than the linear array and 
the test time needed for a full spectrum was quite long. The linear array method was therefore preferable in cases that the 



Wavelength (Angstrom) Wavelength (Angstrom) 

Fig. 6. Measured and calculated Raman Q-b ranch spectra of nitrogen at two different temperatures: 
(a) T=290 K; (b) T=1563 K. 


6 RESULTS 


61 Rayleigh scattering results 

Temperature, velocity, and species number density measurements made with the cw Argon-Ion laser on a different 

rocket engine were reported in an earlier paper. 4 7 Preliminary results of pulsed velocity measurements on the above 

described engine are shown in Fig. 7. Axial and radial velocities were compared against flow field predictions obtained 

using the RPLUS-code. 6 * * * * * * * * * 16 Temperature and number density estimates were not reliable due to uncertainty in the species 

composition in the exit plane of this rocket. The presence of oxygen in the core caused a distortion of the pure Gaussian 

shape based on scattering water molecules that was the basis of the model function. Vertical error bars included both the 

random statistical error associated with photon noise as well as bias errors caused by uncertainty in the scattering angles 

and Fabry-Perot mirror spacing. The horizontal error bars represent the spatial sampling extent of the probe volume which 

was approximately 3 mm in radial direction. The actual uncertainty in the center position of the probe volume was on the 
order of 0.5 mm. 
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The RPLUS code has the capability of making near-field plume predictions for expansion in a vacuum, yet can not 
accomodate flow expansion into a finite pressure environment. Experimental conditions imposed a finite back pressure on 
the exhaust flow. Since Rayleigh measurements were done parallel to the exit plane at 3 mm downstream, some 
discrepancies were expected at the perimeter of the exhaust jet, but the core was expected to be unaffected. The axial 
velocity profile showed a considerable discrepancy between the measured and predicted values, both at the perimeter and at 
the core. It was believed that the difference in velocity gradient was caused by a stronger mixing between the fuel-film 
cooling and oxygen rich core flow than was modeled. This enhanced mixing would lead to a gradually varying molecular 
weight with radius and accordingly would cause a gradually varying velocity going radially outwards. The measured axial 
velocity in the core was lower than the predicted velocity. These differences were attributed to differences in the species 
composition and energy release in the core. 


The radial velocity profile showed qualitative agreement with the predicted values, especially close to the core. Further 
from the core in radial direction the discrepancy increases. It was expected, however, that the back pressure caused the 
expansion to be confined, whereas the predicted expansion was free, allowing uninhibited expansion to occur. This effect 
was most significant at the outermost point. In general, the radial velocity was well predicted, indicating that the nozzle 
divergence was adequately modeled. This conclusion was also supported by previous measurements obtained using cw 
Rayleigh scattering in the plume of a different nozzle, 7 where similar good agreement in radial measurements was obtained. 


(a) 


i ) 

o 

. j 
Ld 


5000 r 
4000 - 
3000 j- 
2000 - 
1000 - 
0 

iooo 
- 10 



Fig. 7. Comparison of Rayleigh velocity data with numerical RPLUS-code prediction: (a) axial velocity component, 
(b) radial velocity component. 

6 2 Spontaneous Raman Scattering Results 


Spontaneous Raman scattering was used to measure species densities and temperatures. However number densities 
predicted in the exit plane of a full nozzle were on the order of 10 16 cm* 3 , and tests showed that Raman scattering was too 
weak to obtain reliable results over the limited test periods. The nozzle was therefore shortened to a 1.86:1 ratio, and 
Raman measurements were completed in a plane 10 mm downstream of this exit plane. Since the nozzle was severely 
undercxpandcd, it was expected that the back pressure would not have a significant effect. Even at this location, a relatively 
low number density was expected, and a PMT with scanned spectrometer was used. Of all the major species expected 
(hydrogen, oxygen, and water), only oxygen and water vapor generated a detectable Raman signal. The oxygen signal was 
strong enough to allow extraction of accurate data. The water vapor spectrum, due to its spread as the result of high 
temperature and its many Raman lines, could not be clearly defined over the limited (30 sec.) test time. Temperature and 
oxygen number densities were extracted from the Stokes line, with calibration factors obtained from oxygen in air. 
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Measured temperature and oxygen number density profiles are shown in Fig. 8a and b, respectively, at 10 mm 
downstream of the shortened nozzle, along with profiles predicted by the RPLUS code. Species detection at this location 
have implications for the the general quality of combustion chamber mixing and combustion efficiency . Both predicted and 
measured temperature profiles showed qualitative agreement, but the magnitude of the measured temperature was about 
10% higher and increased towards the perimeter of the plume. The local maximum in the predicted density profile at about 
5 mm radial position, was caused by a weak compression wave. Such a wave would not be detectable with the Raman 
scattering diagnostics. The corresponding local maximum in the measured temperature profile was thought to be caused by 
the shear layer mixing and combustion between the fuel-film and the oxygen rich core. This caused a higher temperature 
away from the core. This was also predicted in the model, as was shown in the temperature profile at the exit plane (dashed 
line), where the effects of the compression wave and shear layer coincided at about 7 mm radial location. However, in the 
prediction where expansion in vacuum was assumed, this high temperature peak propagated away from the core and was at 
about 12.5 mm radial location when it reached the plane where measurements were taken. As previously mentioned, the 
presence of the back pressure and the diffuser during the measurements limited the expansion. Vibrational non-equilibrium 
of the exhaust gases at the measurement location could be the cause of the difference between the measured and predicted 
temperatures. The temperature was not measured beyond 7.5 mm radially because no oxygen was detected at those 
locations. 
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Fig. 8. Raman data in the exit plane of a low area ratio nozzle: (a) temperature, (b) 02 number density. 


The measured oxygen number density profile was significantly different than the predicted profile in both magnitude 
and shape. Again, the solid line represents predictions at the 10 mm downstream plane where measurements were taken 
and the dashed line represents exit plane conditions. Both the measured radial gradient and number density magnitude lie 
in between the two predicted locations. It is believed that the difference was caused by inadequacies in injector modeling. 
The model assumed plug flow for injection, with homogeneous concentrations and temperature throughout, and with 
constant velocity. Mixing between this core flow and the fuel film takes place in the shear layer originating behind the bluff 
fuel-film insert. It was concluded from the measurements, however, that the injector flow was stratified, with combustion at 
the perimeter, where good mixing between oxygen and hydrogen occurred. Measurements also suggest that hydrogen did 
not penetrate the core flow, leaving the final injection flow very much stratified in number density, and to a lesser degree in 
temperature and velocity. Mixing took place between the oxygen rich injector flow and fuel film, but was not sufficient to 
fully mix over the short distance to the throat. The core remained oxygen rich. This was reflected by the measured profile, 
where the number density of oxygen in the core at 10 mm downstream was almost identical to the predicted number density 
at the exit plane. The same reason caused the measured radial profile to be much steeper, further magnified by the fact that 
the experimental radial expansion was slower than the predicted expansion in vacuum, as was mentioned earlier. 
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A check of the experiment with mass conservation was not possible since water number densities could not be measured, 
even though water did show up as a very weak spectrum. In agreement with predictions, hydrogen could not be detected due 
to the low number densities. Because of its low molecular weight and the fast expansion at the perimeter of the plume in 
vacuum, number densities at the measurement location were calculated to be of the order of 1016 cm-3. This was below the 
detectable limit. 


7 CONCLUDING REMARKS 

Rayleigh and Raman scattering have successfully been applied to determine velocity, temperature and species number 
densities in the plume of a low thrust hydrogen-oxygen rocket installed in a high altitude chamber. Both cw and pulsed 
Rayleigh data and pulsed Raman data were obtained. Much of the transmission of incident and collected radiation was 
facilitated by the use of optical fibers, which were an integral part of the diagnostics equipment. They permitted a reduction 
of stray and scattered light by reducing the number of optical windows. These fibers also provided flexibility in design and 
installation. Optimizing the fibers for diameter and radiation wavelength provided a powerful and multifaceted diagnostics 
facility. 

The measured radial velocity profile was closely predicted by a viscous nozzle code, RPLUS. The measured axial profile 
significantly differed from predictions, both in magnitude and distribution. The cause of this discrepancy was attributed to 
inadequate modeling of the turbulent shear layer mixing in the combustion chamber. Temperature and number densities 
were extracted from the Rayleigh measurements based on assumptions about the species composition. Raman scattering 
data showed that the oxygen number density was different from values predicted using the CFD code. Because oxygen has a 
strong influence on the shape of the Doppler signal, these Rayleigh temperature and number density results for the rocket 
under investigation were not reliable. 

The measured temperature profile determined from Raman scattering spectra was approximately 10% higher than 
predicted. This may be the result of vibrational non-equilibrium. Also the expansion was modeled to occur in a vacuum 
The experimental expansion was much more confined, giving a higher temperature than predicted. Discrepancies also exist 
between the measured and predicted oxygen number density. This was attributed to the initial conditions in the code, which 
included a homogeneous plug flow entering the combustion chamber. In addition, difficulties in accurately modeling the 
turbulence in the shear layer between the fuel-film and the core region was thought to contribute to the differences 

Hydrogen spectra could not be detected and water spectra were too weak to yield accurate results. Raman scattering 
therefore is a more powerful diagnostic technique in a high number density environment, such as in the combustion 
chamber. Vibrational-translational relaxation problems in a supersonic flow biased the temperature measurement toward 
the vibrational temperature. However, with better detection equipment and a longer available test time, vibrational non- 
equilibrium effects in the nozzle flow could be studied. 
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